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> 1 INTRODUCTION 



In a recent study, Adelberger et al (2003, hereafter A03) 
have analyzed a sample composed of 8 high-resolution QSOs 
spectra at 3.1 < z < 4.1 and spectroscopic redshifts for 431 
Lyman-break galaxies (LBGs) at slightly lower z, to study 
the intergalactic medium (IGM) close to high-z galaxies and 
detect imprints of galaxy-IGM interplay. Comparing the lo- 
cation of galaxies to the absorption lines in QSOs spectra, 
the authors found that the mean Lya transmitted flux gen- 
erally increases with the line-of-sight impact parameter be- 
tween the foreground galaxies and the background QSOs 
(Ar), reaching the mean transmissivity of 0.67 at Ar ^ 6 
comoving Mpc; hovever, the trend is reversed within 
about ~ lh~^ comoving Mpc of the galaxy. This result indi- 
cates that the IGM contains less neutral hydrogen close to 
LBGs, and hence is more transparent to Lya photons, than 
at large distance from them. 

Similar observations at lower redshift show an opposite be- 
havior. Lanzetta et al (1995), Chen et al (1998) and Pas- 
carelle et al (2001) find that galaxies, at z < 1, with an 
impact parameter Ar ^ 180/i^^ comoving kpc are much 
more often associated with Lya absorption lines than those 
at larger impact parameters. Furthermore the result of 
VLT/UVES study of the Lya forest in the vicinity of the 
LEG MS1521-CB58 at z = 2.724 (Savaglio, Panagia & 



ABSTRACT 

Recent observations have shown that the intergalactic medium (IGM) is more trans- 
parent to Lya photons close to Lyman Break Galaxies (LBGs) than at large distance 
from them, i.e. a proximity effect. Cosmological simulations including winds from 
LBGs have been so far unable to explain this trend. By coupling such simulations 
with the radiative transfer code CRASH, we investigate whether the addition of the 
ionizing radiation emitted by LBGs can increase the transmissivity by decreasing the 
neutral hydrogen fraction in the inner Mpc of the galaxy halo. The transmissivity as 
a function of distance is roughly reproduced only if LBGs are identified with dwarf 
galaxies (with masses ^ IO^Mq) which are undergoing a vigorous {5QMq yr~^) burst 
of star formation. Similar star formation rates in larger galaxies are not sufficient to 
overwhelm the large recombination rates associated with their denser environment. If 
so, photoionization partly reconciles theory with observations, although we discuss a 
number of uncertainties affecting both approaches. 
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Padovani 2002) shows an absorption excess close to the 
galaxy. 

AOS suggest that the relative lack of neutral hydrogen 
at small distances from LBGs is possibly due to the dy- 
namical feedback from supernova-driven winds, which may 
be able to accelerate the surrounding gas up to velocities 
exceeding the escape velocity of the galaxy potential well, 
driving plausibly the intergalactic material out to a distance 
> 0.5h~^ comoving Mpc. 

This scenario has been investigated with numerical 
simulations by different groups: Croft et al (2002), Br- 
uscoli et al (2003), and KoUmeier et al (2003). Despite 
the completely independent prescriptions adopted to model 
supernova-driven winds in these studies, none of them were 
able to reproduce A03 results on the entire distance range. 
Unless feedback schemes are still far too crude in SPH sim- 
ulations, the discrepancy indicates that the energy injec- 
tion from supernova explosions is not able to displace the 
IGM around galaxies to distances ^ 0.5h~^ comoving Mpc. 
The most likely explanation is that the ram-pressure of the 
accreting/infalling gas slows down and effectively confines 
winds to relatively modest distances from the explosion sites 
(BruscoU et al. 2003). 

Local photoionization by galaxies themselves provide 
the only known alternative mechanism able to decrease the 
neutral hydrogen fraction in the LBG environment. The first 
attempts to account for local photo-ionization in numerical 
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simulations have been performed using simplified analytical 
techniques and/or ad-hoc prescriptions (AOS, Kollmeier et 
al. 2002, Croft et al. 2002); these studies have lead to con- 
clude that this process cannot account for the observed IGM 
transparency. However, fully sclf-consistcnt and detailed ra- 
diative transfer calculations are required before a final con- 
clusion is drawn. 

In this Letter we investigate the role of local pho- 
toionization using high resolution numerical simulations per- 
formed with the radiative transfer code CRASH ; the full de- 
scription of the code and of the tests made are given in 
MaselU, Ferrara & Ciardi (2003). 



2 SIMULATIONS 

We have performed radiative transfer simulations in order 
to evaluate self-consistently the proximity effect due to the 
ionizing radiation emitted by high redshift galaxies on the 
surrounding IGM. Wc use the numerical code CRASH (Maselli 
et al. 2003) which allows to calculate the effects of ioniz- 
ing radiation on a precomputed cosmological density field, 
specifically obtained from SPH simulations fully described 
elsewhere (Marri et al. 2003) and already used by Bruscoli 
et al. (2003). In the following we give the details of the two 
simulations. 

2.1 Hydrodynamics 

The simulation has been carried out with an improved ver- 
sion of the SPH code GADGET2 (Springcl et al. 2001), de- 
signed for the detailed study of galaxy formation (Marri & 
White 2002). Such multi-phase (MSPH) scheme is based on 
an explicit separation of protogalactic gas into diffuse and 
dense (star forming) components and includes new prescrip- 
tions for a realistic treatment of stellar feedbacks. With re- 
spect to the standard GADGET2 version, this new imple- 
mentation is able to suppress star formation, to reheat cold 
gas and to drive outflows from galactic disks. The presence 
of a UV background (UVB) produced by QSOs and filtered 
through the IGM is also included. The shape and amplitude 
of the UVB are taken from Haardt & Madau (1996) and its 
impact on the IGM temperature and ionization structure is 
calculated as in Katz et al. (1996). 

The simulation uses 128^ particles in a W.5h~^ comov- 
ing Mpc box and assumes a ACDM cosmological model with 
Qo = 0.3, Qa = 0.7, 06 = 0.04 and h = 0.7; the initial 
density perturbations power spectrum is cluster-normalized 
(us = 0.9) and periodic boundary conditions arc assumed. 
Galaxies arc identified by applying the HOP group finding 
algorithm (Eiscnstcin & Hut 1998) to the baryonic compo- 
nents of the simulation. Each galaxy is assigned a star for- 
mation rate (SFR) calculated by considering the properties 
of all particles within the group. For each of them wc first 
check if the conditions under which a gas particle can form 
stars are satisfied and we calculate the rate at which its gas 
is converted into stars. For more details we defer the reader 
to Marri & White (2002). 

For consistency sake, wc use the same output as in 
Bruscoli et al. (2003). There we showed that, although on 
large scales the simulation reproduces a number of statisti- 
cal properties of the IGM, it does not match the transmitted 



Lya flux inside 1 h^^ comoving Mpc as measured by A03. 
Such an output has now been post-processed using CRASH to 
self-consistently include the effects of local photoionization 
produced by the galaxies identified in the computational vol- 
ume. 

2.2 Radiative Transfer 

The radiative tranfer (RT) simulation has been performed 
on the outputs at z = 3.27 of the MSPH simulation above, to 
allow a direct comparison to the observational data available 
around this redshift ( A03) . 

The RT is simulated running the 3-D RT code CRASH , 
whose description is fully reported in Maselli et al. (2003). 
The code works on precomputed H/He density fields as- 
signed in a 3-D cartesian grid, and calculates the time evo- 
lution of the temperature and of the ionization structure. 

To include the effect of local photoionization on the 
cosmological configuration obtained from the MSPH simu- 
lation, it is necessary to map the particle data on a 3-D 
cartesian grid. We choose a grid composed of 128^ cells, 
which corresponds to a spatial resolution of Ax ^ 82 kpc 

comoving. At each grid point we assign a value for 
the mass density (p), temperature (T), neutral hydrogen 
fraction (a^no = nHo/nn) and free electrons fraction (xe = 
Mc/(^^H + ^^Hc)), performing a standard SPH smoothing us- 
ing the 32 neighbour SPH particles closer to the grid point. 
Wc next locate the 398 galaxies identified inside the simu- 
lation box on the cartesian grid, converting the lagrangian 
spatial coordinates into cartesian coordinates. 
The emission properties of the galaxies have been derived 
using the Starburst99 code (Leitherer et al. 1999), a web 
based software and data package designed to model spec- 
trophotometric and related properties of star-forming galax- 
ies. We have assumed a continuous star formation law and 
a Salpeter IMF with a — 2.35, normalized in a mass range 
[1,100]Mq; we also take Zq/5 as a characteristic metallicity 
for all the sample. These values are suggested by the analy- 
sis of the LBGs UV continuum and the H^ and H^j emission 
lines by Pettini et al. (2001). As the physical connection 
of LBGs with galaxy type/mass is unclear (Primack et al. 
2003), we assume for all gala^xies in the box similar scalings. 
Running Starburst99 under the above conditions, we obtain 
the galaxy SED corresponding to the derived SFR of each 
galaxy. We finally assume a value for the ionizing photon es- 
cape fraction of 50%. This value represents the upper limit 
of the allowed range (see discussion in Ciardi, Bianchi & 
Ferrara 2002) from a set of observational data. 

The computed grid data are assigned as initial condi- 
tions for the CRASH run. The number densities of hydrogen 
and helium nuclei are calculated assuming mass fractional 
abundances /if=0.76 and /ife=0.24, respectively. Tempera^ 
ture and hydrogen ionization fractions are initialized as in 
the MSPH simulation output; for helium we assume ioniza- 
tion fractions = and a;jjo++ = (s^e— a;Ho+~a;jj+)/2. 
The physical duration of a simulation run has been set 
to = 6 X 10^ yrs, i.e. several times the estimated box- 
crossing-time for the ionization front produced by a typical 
source in the volume. 

Nearby galaxies the mean gas temperature is above 10^ 
K, due to the energy injection into the IGM produced by 
outflows. The cooling time of such gas in the simulation is 
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Figure 1. Neutral hydrogen fraction maps in a slice of about 4h~^ comoving Mpc by side, centered on the most massive and most 
luminous galaxy in the computational volume. The three slices represent, in a logarithmic scale, the final outputs of run A (left panel), 
run B with the galaxy SFR derived from the MSPH simulation (central panel), and run B with the SFR boosted by a factor of 10 (right 
panel). 



longer than t^, with a minimum value equal to 1.5x 10* yrs. 
As the temperature does not change significantly during the 
RT simulations, its evolution has not been solved for to re- 
duce the computational time. 

Two different RT simulations have been run. In the first 
one only the UVB radiation is included (run A), while in 
the second one, the radiation from all galaxies in the sim- 
ulation volume is added to the UVB (run B). In this way 
it is possible to compare consistenly the two simulations to 
infer the role of galaxies and distinguish the contribution to 
the proximity effect due to local photoionization. As consid- 
erable uncertainty is present in the determination of SFRs 
of LBGs, we also studied a case similar to run B in which 
however SFRs are boosted ad hoc to match the observed 
values. 

The total ionizing energy, injected in the computational 
volume, is discretized in Np — 2 x 10** photon packets for 
run A and in Np = 6 x 10* packets for run B: roughly, this 
corresponds to 10® photon packets per galaxy and 2 x 10* 
photon packets for the UVB. 



3 RESULTS 

As an example, Fig.l shows a comparison among the final 
outputs of run A (left panel) and run B (central and right 
panels) . The slices represent the logarithm of the neutral hy- 
drogen fraction on a slice of a zoomed region of about 4:h~^ 
comoving Mpc by side, centered on the most massive galaxy 
in the computational volume. The two maps from run B 
have been derived by assigning to such a galaxy different 
SFR: in the central panel the SFR is about 30 Mq yr~^ 
as derived from the MSPH simulation, while in the right 
panel it has been boosted by a factor of 10. The effect of 
the local photoionization produced by the galaxy is clearly 
visible in both cases and depends strongly on its luminosity, 
as expected. In the central panel the imprint of the galaxy 
on the surrounding gas is feeble and the content of neutral 



hydrogen in a surrounding region of 0.5 h~^ comoving Mpc 
by side, is reduced on average by factor of 2. The luminosity 
associated with a SFR of 30 M0 yr~^ is too low to affect 
significantly the high density environment of such a massive 
galaxy (M =1.5 x IO^Mq/i-^). When the SFR is boosted by 
a factor of 10 (right panel), the region affected by the galaxy 
UV flux extends up to a distance of about lh~^ comoving 
Mpc. In the inner region the neutral hydrogen fraction is 
depleted at least by an order of magnitude. In the outer re- 
gions the effect of the local photoionization is progressively 
reduced and vanishes at distances lh~^ Mpc, where the 
UVB dominates the local UV flux. The decrease in the neu- 
tral hydrogen fraction visible in other regions far from the 
center is produced by other galaxies, included in run B as 
ionizing sources. 

To estimate the radius of the "sphere of influence" pro- 
duced by the local emission from a typical galaxy we study 
the distribution of the ratio between the UV flux produced 
by the galaxies and the UVB flux, TZ — Fgai/Fbkg, in the 
box. Fig. 2 shows the histogram derived from the final output 
of run B. The distribution has a mean value TZ — 2.64 x 10~^ 
and a median value TZm = 1.83 x 10~^. These values are to 
be considered as lower limits to the contribution of galaxies 
to the UVB, as larger boxes would be required to prop- 
erly estimate such quantity. The volume in which TZ > \, 
i. e. where the local emission by galaxies dominates the fore- 
ground UVB, is about 4.8% of the total computational vol- 
ume. This fraction corresponds to a mean radius of infiuence 
for a typical galaxy Ri, given by: 



4 iV,(7^ > 1) x K 



37r 



1/3 



0.4/i"^Mpc 



(1) 



where is the number of cells where 7?, > 1, 14 is the vol- 
ume of a single cell and Ng — 398 is the number of galaxies 
in the simulation box. 

To allow a direct comparison of the simulation results 
with the observational data we construct synthetic Lya 
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Figure 2. Distribution of the ratio between the UV flux produced 
by the galaxies and the UVB flux, TZ = Fgai/Fikg, in the box, as 
derived from the final output of run B. 



spectra as described in Bruscoli et al. , by tracing ran- 
dom LOS through the box. We then estimate the mean Lya 
transmitted flux as a function of the impact parameter Ar, 
computing the average on all the pixels at a distance from 
a galaxy in a given interval centered on Ar. The results are 
plotted in Fig. 3 and compared with the data, represented 
by the black points. 

The three panels show the mean Lya flux computed on 
different samples of galaxies: all the 398 galaxies (left panel), 
galaxies with mass above 2 x 10^° Mq (central panel) and 
those with mass below 9.3 X 10* Mq (right panel). 

Solid, dotted and dotted-dashcd lines arc derived from 
the outputs of run A, run B and from a simulation analogous 
to run B, but with SFRs boosted ad hoc. The difference 
between run A and run B in the left panel is marginal: we 
conclude that the mean effect of the local photoionization 
on the Lya transmitted flux is negligible when all gala^xies 
are included in the analysis. 

Massive galaxies are the best candidates for LBGs due 
to their high luminosity and clustering properties. The SFR 
derived for the most massive galaxies in the MSPH simu- 
lation are in the range 10 30 Moyr"'^, the highest values 
in the simulation. Due to their high luminosity, one would 
expect them to produce a strong impact on the ionization of 
the surrounding gas. However this is not the case. Looking at 
the solid line in the central panel, obtained by neglecting the 
local emission, we note that the high density characterizing 
the environment of massive galaxies suppresses the mean 
transmitted flux with respect to the mean trend obtained 
for all the galaxies in the simulation (solid line in the left 
panel). The UV radiation emitted by these galaxies, using 
the SFR derived from the MSPH simulation (dotted line in 
central panel) , is not strong enough to enhance significantly 
the transparency of the surrounding IGM, because of the 
high recombination rate in the denser environment. Higher 
values of SFR are necessary to produce a significant increase 
of the transmitted fiux nearby massive galaxies. This can 



be seen from the boosted case, i.e. with SFR in the range 
100 300 Mq yr~^. Although on the high side, the values 
are still plausible according to observations (Shapley et al. , 
2001). Despite the sharp increase in the mean Lya transmit- 
ted flux, yet this simulation does not match the A03 data. 
This result suggests that the mean transmissivity strongly 
depends on the galactic enviroment. The right panel shows 
the trend for the Lya transmitted flux obtained selecting 
galaxies with mass ^ 9.3 x 10* Mq. In this case the trans- 
missivity of the gas at Ar < Q.Th^^ Mpc has an opposite 
trend with respect to the mean (left panel) , and seems to fol- 
low the one observed by A03. The increasing fiux at smaller 
distance from the galaxy reflects the lower density in these 
regions. Here the UVB flux can counteract the recombina- 
tions which occur on longer time scales. Nevertheless the 
SFRs associated with these galaxies are too small (of order 
O.IMq yr~^) to further increase the transmitted flux. Boost- 
ing the SFR up to values around 50 Mq yr~^ (dashed-dotted 
line), the increase in {F) is signiflcant and closely approches 
the observed value in the innermost region. By no means 
though the flux around 1 Mpc can be matched. 



4 CONCLUSIONS 

A recent study on the statistics of the Lya forest in the 
vicinity of foreground galaxies has revealed an unexpected 
relative lack of neutral hydrogen in the inner Mpc co- 
moving from LBGs at « ~ 3. High energy supernovardriven 
winds able to displace the surrounding gas out to distances 
greater than 0.5/i~^ have been proposed as a viable expla- 
nation. However, several numerical studies failed in repro- 
ducing the observed effect of such superwinds via (M)SPH 
cosmological simulations. According to the simulations, the 
velocity of the gas closest to galaxies varies around a char- 
acterisitc value of 40 km s~^ for small objects, and around 
150 km s~^ for the biggest ones. These values seem to un- 
derestimate the high velocity gas, at 600 km s~^, measured 
by A03 at a scale of a few kpc from the galaxy center. Nev- 
ertheless, the possibility for the wind to break through the 
galactic halo, and into the IGM up to distances greater than 
lh~^ comoving Mpc (even at such high velocity) seems quite 
unphysical if one takes into account the kinetic energy used 
to counteract the pressure of the accretion flow. 

Motivated by these reasons, we have studied the role 
of local photoionization in determining the Lya absorbers 
statistics in the vicinity of galaxies, via numerical radiative 
transfer simulations. We post-processed an output aX z = 
3.27 of a multiphase SPH simulation, in which a consistent 
treatment of supernova feedback is implemented, running 
the code CRASH to account for the impact of the UV flux 
emitted by galaxies on the simulated IGM. We have derived 
the synthetic mean Lya transmitted flux, as a function of 
the distance from a typical galaxy, for different samples of 
galaxies. We found that on average local photoionization has 
a negligible impact, i.e. when all the galaxies present in the 
box are considered in the flux distribution derivation. 

Local phoionization can affect significantly the trans- 
parency of the surrounding IGM only for particular condi- 
tions of the galactic environment. We find that in order to 
reproduce the trend inferred by AOS, the gas around LBGs 
should be less dense than that surrounding massive galax- 
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Figure 3. Mean Lya flux at diflferent impact parameters Ar, computed as the average values among all the pixels at a distance from 
a galaxy in a given interval centered on Ar. The curves in the three panels are derived using different samples of galaxies: all the 398 
galaxies (left panel), galaxies with mass above 2 X lO^'' Mq (central panel) and those with mass below 9.3 X 10* Mq (right panel). 



ies seen in the simulation. The transmissivity as a function 
of distance is roughly reproduced only if LBGs are identi- 
fied with dwarf galaxies (with masses < IO^A/q) which are 
undergoing a vigorous (SOM© yr~^) burst of star formation. 
Similar and higher star formation rates in larger galaxies are 
not sufficient to overwhelm the large recombination rates 
that are occurring in their denser environment. This con- 
clusion is somewhat at odd with the common wisdom that 
LBGs are predominantly massive systems, although differ- 
ent and equally valid arguments tend to identify them with 
the bursting population of dwarfs that we find necessary to 
explain the observed gas transparency. LBGs can be inter- 
preted as outbursting dwarf galaxies (Bullock et al. 1999) if 
their abundance is due to an increased collision rate at high 
redshift; although many of the simulated collisions have rela- 
tively small masses ( ^ 10^*^1/0), they tend to cluster about 
large-mass halos. They therefore exhibit strong clustering, 
similar to that observed. 

The solve the puzzle, more data are needed. The poor 
statistics of AOS sample at small impact parameters, errors 
in measured redshifts and galaxy positions make the derived 
transmitted flux quite uncertain. 
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